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Based on Run I data we present a comprehensive analysis of Higgs couplings. For the first time 
this SFitter analysis includes independent tests of the Higgs-gluon and top Yukawa couplings, Higgs 
decays to invisible particles, and off-shell Higgs measurements. The observed Higgs boson is fully 
consistent with the Standard Model, both in terms of coupling modifications and effective field the¬ 
ory. Based only on Higgs total rates the results using both approaches are essentially equivalent, 
with the exception of strong correlations in the parameter space induced by effective operators. 

These correlations can be controlled through additional experimental input, namely kinematic dis¬ 
tributions. Including kinematic distributions the typical Run I reach for weakly interacting new 
physics now reaches 300 to 500 GeV. 
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I. INTRODUCTION 

The discovery of a light narrow Higgs boson [T] in 2012 [2] was a triumph of particle physics. While its statistical 
significance was largely driven by the peak in the 77 invariant mass spectrum, already Run I of the LHC allowed 
ATLAS and CMS to perform a large number of tests of the nature of the observed resonance. One of these tests is 
the analysis of Higgs couplings relative to their Standard Model values. No significant deviations from the Standard 
Model properties were observed in the Higgs production and decay rates. However, we need to keep in mind that 
these constraints are at a numerical level where in typical weakly interacting models for new physics |3) we would not 
expect significant deviations, either. 

One of the main physics programs defining the upcoming LHC runs will be a comprehensive precision analysis of 
Higgs properties. It will eventually utilize up to 3000 fb _1 of integrated luminosity, covering a wide range of production 
and decay channels and observables. Technically, the Higgs coupling analysis from Run I is based on comparably 
simple total cross sections and branching ratios; this simple structure of underlying measurements allows us to limit 
the interpretation to an independent variation of all Higgs couplings in the Standard Model without missing much of 
the experimental information. 

Planning for the corresponding Run II analysis we need to implement major extensions to the Higgs coupling 
analysis. The reason is that the expected wealth of measurements, including kinematic distributions, will probe 
modifications of the Standard Model Lagrangian which are not captured by a simple shift in Higgs couplings.* 
Including kinematic distributions in the analysis of the Higgs couplings suggests to describe the theoretical basis in 
terms of an effective field theory. It also implies that the interface between experiment and theory has to be re-defined. 
Eventually, at least three questions need to be addressed by the Higgs measurements: 

1. do any of the observed Higgs couplings show deviations from their Standard Model predictions? 

2. can possible deviations be linked to a mass scale through an effective theory of the electroweak Lagrangian? 

3. what does the Higgs sector tell us about specific new physics models? 

These questions are closely tied to each other, but the complex structure of correlations between the experimental 
measurements and between theoretical predictions render a single answer to all three questions unrealistic. In the 
past, SFitter has been one of the driving forces to answer the first SHU and third [3] questions. The second question 
has been tackled by Higgs specialized tools, for example presented in Refs. ei m - A shift from the analysis of 
Higgs couplings based on total rates to an effective field theory of the Higgs sector is the appropriate step to include 
kinematic distributions in the Higgs fit. 

In this paper we first update the Higgs coupling analysis to include the full Run I data set. One focus of SFitter 
is the appropriate treatment of theoretical uncertainties, which will also play a major role in this paper. Second, 
we present a first SFitter study of the electroweak effective held theory, which will in the future allow us to treat 
the three above questions in the same framework. With that purpose we include in the analysis several differential 
distributions and we study their implementation and impact. The Run I constraints in terms of the accessible new 
physics scale stay below the TeV range, which is consistent with a 20% precision on coupling analyses and supports the 
statement that within known theoretical frameworks we do not expect to see deviations in the Higgs couplings during 
Run I. Finally, we include off-shell measurements of Higgs production HU in the SFitter analysis and estimate 
their current impact. While it has been established that off-shell measurements will not serve as model-independent 
measurements of the Higgs width m, we can establish their role in hypothesis-based Higgs studies. 


This does not mean that the Higgs coupling analysis should be abandoned entirely. For reasons discussed in this paper we will keep it 
in the SFitter framework to search for physics beyond the Standard Model which merely shifts Higgs couplings. 
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A. Experimental input 


In all SFitter Higgs coupling analyses the experimental input is not the published set of Higgs signal strengths, 
but the number of signal and background events for each analysis. This allows us to independently study statistical, 
systematic, and theoretical uncertainties and include them in a profile likelihood analysis largely independent of the 
ATLAS and CMS assumptions. As a matter of fact, this independence is our motivation to maintain the SFitter 
effort in spite of more and more advanced experimental Higgs coupling analyses. 


As experimental input we use the following Higgs searches and measurements as published by the ATLAS and CMS 
collaborations: 


production/decay mode 

ATLAS 

CMS 

H —► WW 

Ref. [13] 

Ref. [13 

H -s- ZZ 

Ref. [15] 

Ref. [16] 

H —i > 77 

Ref. [13 

Ref. [ig 

H —> rf 

Ref. [TO] 

Ref. [HI 

H^-bb 

Ref. [23 

Ref. [22] 

H ^ Z-i 

Ref. [23] 

Ref. [13 

H —>■ invisible 

Ref. [MU] 

Ref. mm 

ttH production 

Ref. [HI S3 Ref. [131 [32] [33] 

kinematic distributions 

Ref. E3[34] 

off-shell rate 

Ref. [35] 

Ref. |36| 


From all these analyses we extract the number of observed, signal and background events after appropriate cuts. The 
several categorizations in the experimental searches listed in the above table lead to the 159 measurements that we 
include in the rate based analyses. In Sec. |III D] we add 14 extra measurements from kinematic distributions. Finally, 
the off-shell distributions considered in Sec. IIVI contribute with 37 additional measurements. 


We will show Higgs coupling analyses in the SFitter framework starting with one universal modification, moving 
to five tree level SM-like couplings, and then allowing for additional contributions to the Higgs-photon and Higgs-gluon 
loop-induced couplings. A new channel which we did not include in previous Higgs coupling analyses is the direct 
measurement of a modified top-quark Yukawa coupling in ttH production. This is necessary to identify new physics 
contributing to the effective Higgs-gluon coupling. While there is not yet enough sensitivity to properly observe the ttH 
production channel with SM coupling strength, the current searches do provide upper bounds on the production and 
decay rates. Both experiments have looked for this channel as part of their di-photon final-state analysis [TTt lTBI [31] : 
moreover, CMS has published additional searches based on multi-lepton final states arising from Higgs decays into 
WW*, ZZ* and rf [32] as well as a dedicated bb analysis [33]. These ttH measurements allow us to separate an extra 
new contribution to the Higgs-gluon coupling, leading to the analysis with seven independent coupling modifications. 
After setting the two Higgs-weak-boson coupling modifications equal, these seven modifications correspond to the 
relevant parameters in the non-linear effective Lagrangian expansion [37], if we restrict them to the leading terms 
following Ref. 38]. In the last step we will also include Higgs searches to invisible particles, i.e. generating missing 
transverse momentum. This will allow us to add an extra modification to the analysis, accounting for possible Higgs 
invisible decays. All coupling analyses in the traditional SFitter framework rely on the same measurements listed 
above. 


For the dimension-6 operators, which not only change the coupling strengths but also the momentum dependence 
of the vertices leading to modifications of kinematic distributions, we make use of several differential distributions. 
These are the transverse momentum distributions of the gauge boson in VH, H —> bb production [21] for all 0, 1 
and 2-lepton final states, and the A <j>jj distribution in H +2 jets with di-photon decays of the Higgs boson (3J. In 
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Sec. IIII Cl we will discuss in detail how we include the crucial information from kinematic distributions in addition to 
total rates. 

Finally, we will for the first time include off-shell Higgs production rates [35], [SB] in a global fit of Higgs properties. 
We will analyze their effect to the usual coupling determination and quantify their potential to constrain the Higgs 
total decay width. 


B. Fit setup 


The most relevant shortcoming of Higgs production at the LHC is that there is no method which allows us to 
directly access the Higgs width if we only consider on-shell data. For the traditional Higgs couplings analysis this 
means that we will be sensitive to all observed Higgs decays, including decays to invisible particles passing through 
the detector without leaving a trace. Because the total width entering the coupling measurement is a sum of partial 
widths, we can easily construct a lower limit to the total width by summing the partial widths induced by the observed 
Higgs couplings. 

Once we add off-shell Higgs data we will start being sensitive for example to Higgs decays to hadrons or more than 
two states. The reason is that off-shell Higgs rate measurements add a possible upper limit on the Higgs width, as 
discussed in detail in Sec. m This limit is currently weak, so for the present section we identify the sum of observed 
partial widths with the total width, as usually done in SFitter. The one additional assumption we need to make 
is generation universality, meaning that the relation between the mass and the Yukawa coupling for up-type quarks, 
down-type quarks, and charged leptons is universal for the second and third generations each. For example, the charm 
Yukawa coupling is shifted by g c = gf M (l + A t ) in the usual SFitter conventions which will be reviewed in Eq.(2.11. 

With these assumptions we know what would happen if there was a sizeable unobserved contribution to the total 
width: in the ideal case of all measurements otherwise in agreement with the Standard Model and equally constraining, 
all measured Higgs couplings based on the underestimate of the total width will appear to be too small by a universal 
factor. The experimental challenge would be to distinguish such a scenario for example from a Higgs portal where 
the SM Higgs state mixes with a heavy additional scalar f3S] ■ In this situation the off-shell Higgs rate measurement 
will be extremely useful. 


For all experimental uncertainties we assume a Poisson distribution for the statistical error of the rate measurements 
and Gaussian distributions of associated nuisance parameters. For the systematic uncertainty distribution this is 
justified as long as it is resolved with help of some kind of measurement outside the analysis we are considering. For 
theoretical uncertainties a Gaussian modeling might be technically convenient, but it needs to be established that it 
is at least conservative m■ We will give a detailed discussion of this issue in Sec. 

The usual determination of model parameters in SFitter proceeds in two steps. First, we construct Markov chains 
which probe the multi-dimensional space of Higgs couplings. These chains cover typical values of A x £ [—5,7]. The 
error treatment includes all experimental and theoretical error bars with full correlations. From these chains we for 
example compute the 2-dimensional correlations shown later in this work. In a second step we can focus on the 
SM-like solution, defined as A x > — 1 for all couplings, and compute the 68% and 95% confidence-level (CL) interval 
of the log-likelihood for each individual Higgs coupling. 

In the Higgs analysis presented in this paper we already derive 1-dimensional profile likelihoods and the correspond¬ 
ing 68% CL interval from the Markov chains. Unlike for the 2-dimensional correlations we limit our 68% CL analysis 
in the individual couplings to the SM-like solution. The interval is defined as the range of coupling modifications 
or Wilson coefficients covering 68% of the integrated profile log-likelihood. Because in particular for non-Gaussian 
distributions this definition allows for a simultaneous shift of both limits and is hence not uniquely defined, we also 
require the value of the profile log-likelihood to coincide on both sides of the error band. 


IIC 
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II. HIGGS COUPLINGS 


One question we can ask in the Higgs sector is: how well does the Standard Model describe all available Higgs data 
at the LHC? There are (at least) two ways of answering this question: first, we can compute an over-all confidence 
level of the Standard Model given all available Higgs data, possibly combined with electroweak precision measurements 
etc. The answer to this question is statistically well defined, but to give us useful information we need to carefully 
analyze the pulls of different measurements. 

Second, we can measure the parameters in the renormalizable Higgs Lagrangian: we start from the quadratic and 
quartic terms in the Higgs potential, which can be exchanged for the Higgs mass and the vacuum expectation value 
v. This relation can eventually be tested in measurements of the Higgs self-coupling. In addition, we can measure 
the Higgs coupling to each Standard Model particle. The Lagrangian underlying this measurement consists of the 
Standard Model operators with free couplings. These couplings can be extracted from LHC rate measurements, 

9x = gl M (1 + A x ) 

g 7 = g* u (1 + Af + A 7 ) = 57 sm (1 + Af+ NP ) 

9g = 9g M (1 + A® M + A g ) = g s g M (1 + A® m+np ) , (2.1) 

where the modifications of the tree-level couplings appearing in the Standard Model loops are encoded into A®^ 
while extra possible new physics contributions are included in A 7iS . In general, the loop-induced couplings have a 
non-trivial momentum dependence; for the Higgs couplings measurement we assume that all three external momenta 
are fixed, reducing the coupling to a single number. In terms of a Lagrangian we can write our hypothesis as 

£ = £sm + A w ginwH W^Wg, + A^ nizH ^ A/ —-H (Jr/l + h.c.) 

Zc w V 

T,0,t 

TT TT 

+ A q Fq — GnjyG^ 1 ' + A~Fa — A uiy A^ + invisible decays , ( 2 . 2 ) 

v v 

with c w denoting the cosine of the weak mixing angle. A possible Higgs decay to invisible states could be described 

by a wide variety of Lagrangian terms, but as long as we only search for an invisible branching ratio there is no 

need for further specifications. The contributions to the higher-dimensional Higgs-gluon and Higgs-photon couplings 
are normalized to their Standard Model values Fq and Fa- We use their values for finite loop masses, while their 
normalization is illustrated by the limit F^ x, ' > —> a s /(127r) for a heavy top mass. The invisible decay width can for 
example be generated in a Higgs portal model [39] and will be quoted in terms of an invisible branching ratio. 


The form (1 + A) of the coupling deviations suggests that we will focus on scenarios not too different from the 
Standard Model. Large deviations from the Standard Model should be taken with a grain of salt, because in such 
a situation it is not clear if the kinematic distributions of the ‘Higgs’ signal and the associated detector efficiencies 
are well controlled. Values around A = —2 indicate a switch in the sign of the coupling, which should be checked 
individually. In some figures we show the full range of A values for illustration, but we will limit our interpretation 
to small values of |A| < 0.5. 


The coupling measurement according to Eqs.(2.11 and (2.2) is clearly well defined in the sense that it corresponds to 


a weak scale Lagrangian which we compare to experimental results. Higher-order QCD corrections can be included, 
because renormalizability with respect to the strong coupling is not affected by changes in the Higgs couplings. 
Electroweak corrections cannot be computed in a model with free Higgs couplings, but their impact on LHC rates 
can safely be neglected for Run I data. 


In the spirit of an effective field theory the free couplings ansatz of Eq.(2.1) can be linked to extended Higgs 


sectors [5] : in such models the Higgs coupling measurement will search for modifications to the couplings of the 
SM-like Higgs boson. Using a well defined ultraviolet completion, possible shifts can also be interpreted for example 
using the full set of model parameters in an aligned two-Higgs-doublet model. If we include A t as well as A g in the 
coupling analysis we need to supplement the extended Higgs sector for example with an additional strongly interacting 
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Figure 1: 68% CL error bars on the deviations from all Standard Model couplings of the observed Higgs boson. In this fit 
we do not allow for new particles in the effective Higgs couplings to photons and gluons, A 7 = 0 = A s . The results labelled 
‘SM exp’ assume central values on the Standard Model expectation, but the current data error bars. 


fermion. In addition, the non-diagonal photon-Z-Higgs is missing, but can be trivially added to this ansatz as we 
will discuss. We have shown that a SFitter analysis of the aligned two-Higgs-doublet model - with coupling 
deviations A computed from the underlying parameters — and the weak scale Higgs coupling analysis indeed give 
identical results [8]. Furthermore, neglecting the invisible decays and setting Aw = A z the Lagrangian in Eq.(2.2) 
corresponds to the non-linear effective Lagrangian m of the Higgs sector, but restricted to the leading terms of the 
expansion defined in Ref. [ 581 . 

Our couplings approach only tracks a deviation in the leading coupling of the Higgs boson to each SM particle. We 
only consider the dimension-4 Lagrangian plus non-decoupling dimension-6 operators coupling the Higgs to photons or 
gluons. Strictly speaking, in the linear effective Lagrangian approach this kind of dimension-6 operators also modify 
the Higgs couplings to electroweak gauge bosons. As long as the experimental analysis is limited to total rates, we 
can expect loop-induced corrections for the HWW coupling to be suppressed with respect to shifts in the tree-level 
coupling g w . 


A shortcoming of the approach is that it does not include effects of the additional new particles on the SM 
Lagrangian outside the Higgs sector. For example, it does not link deviations in Higgs couplings to anomalous triple 
gauge couplings. For a comprehensive analysis of all effects of such new states we have to extend the free Higgs 
couplings to an effective field theory based on a non-linear sigma model in the broken phase of the electroweak 
symmetry, see for instance Ref. m- 

Going beyond the original SFitter ansatz indeed means re-writing and extending the Lagrangian by additional 
operators which couple Standard Model fields to the Higgs boson. These operators can be classified by their dimension. 
For example for the Higgs coupling to W and Z bosons this question has been studied independently of the coupling 
strength mmm- 


A. Standard Model couplings 

In the first step we can fit the five tree level Higgs couplings to all Standard Model particles relevant for the LHC 
observations. The result is shown in Fig. [T] The red bars labelled l SM exp’ show results where we have injected 
a Standard Model Higgs signal on top of the background, i.e. the measured rate in each channel is exactly the SM 
expectation, but leave everything else unchanged. They indicate that the observed errors are slightly smaller than 
expected. This is a universal effect of the theoretical uncertainties which we will discuss in detail in Sec. |II C| 


The simplest model, motivated for example by a Higgs portal [39j or a single form factor from a strongly interacting 
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Figure 2: 68 % CL error bars on the deviations A x from all Standard Model couplings of the observed Higgs boson. For the 
loop-induced couplings we allow new contributions to the H 77 coupling only (A ff = 0, left) and to the H 77 and Hgg couplings 
(right). The results labelled ‘SM exp’ assume central values on the Standard Model expectation, but the current data error 
bars. 


Higgs sector gl], consists of a universal coupling modification Ah- Such a coupling modification is constrained to 
around 3% at 68 % CL, in agreement with the Standard Model or A h = 0. Translated into a mixing angle a from 
the Higgs portal whose preferred range at 68 % CL is 

cos a = 1 + A h e [0.93,1.03] . (2.3) 

The second simplest model is a universal coupling modification for the Higgs interaction with the gauge bosons Ay, 
and one modification for the coupling with the fermions A f. Possible ultraviolet completions are given by models 
with additional Higgs multiplets beyond singlets or doublets, where certain combinations allow us to circumvent the 
otherwise strong limits by electroweak precision data |32] . In this case, Ay is still constrained to around ± 6 % at 
68 % CL, while the fermionic coupling shows a reduced precision of around ±12% at 68 % CL, all consistent with the 
Standard Model. 


An independent variation of five Higgs couplings is also in complete agreement with the Standard Model. Again, 
the actual error bars on the A^ are slightly smaller than what we would expect from exact Standard Model values, an 
effect we will discuss in Sec. IIC The 15% measurement of the top Yukawa coupling is driven by the Higgs couplings 
to photons and gluons under the assumption that no new particles contribute to these loop-induced couplings. The 
measurement of the bottom Yukawa benefits from the normalization of all rates, because in the Standard Model the 
total width is largely driven by the partial decay width H —>• bb. 

In addition to the individual couplings we also show the deviations of ratios of couplings. Such ratios are useful to 
remove systematic and theoretical uncertainties. Indeed, we see that the ratio 


9b 

gw 


ffb SM 

Qw 1 


(l + A b / W ) 


(2.4) 


shows a smaller variation than A b alone. The corresponding positive correlation of A b and Ay/ arises from the total 
width in the denominator of the predicted event numbers. 


Because the Higgs decay H —> 77 has been precisely measured at the LHC we can extend the coupling fit by a new 
physics contribution to this loop-induced coupling. Of course, the variations of the Standard Model couplings A b .t 
and Ay/ are consistently reflected in the full A® M+NP . Following Eq.(2.1) this deviation consists of two terms, the 
parametric shift from the Standard Model loops and an additional shift from new physics. The latter is shown as part 
of the coupling measurement in the left panel of Fig. [2] where we present the results of the six parameter analysis. 























































The only modification with respect to the fit with A 7 = 0 is a slight downward shift of the central value of An/. It 
decreases the contribution from the dominant LF-loop and therefore has to be compensated by a small positive new 
physics contribution A 7 ~ 0.13. This also leads to a very slight increase in the uncertainty on Aw and A t . While 
the new physics contribution A 7 based on all available ATLAS and CMS analyses has a one-sigma preference for an 
additional contribution, the combination A 7 M+NP is in perfect agreement with the Standard Model. The error bars 
for A 7 and A® M+NP have the same size, which means that the interference structure between A w and A t breaks any 
strong correlation with A 7 in this fit. 

Finally, we show in the right panel of Fig. [2]the first SFitter Higgs couplings fit including a new physics contribution 
to the effective Higgs-gluon coupling, A g ^ 0. As argued in Sec. [TT| and suggested by our notation, in this seven 
parameter analysis we focus on the SM-like solutions, i.e. small values of the |A X |. As expected, the increase in the 
error bar of A t is dramatic. The central value of A t increases by one standard deviation of the new measurement, 
while A g resides about one standard deviation below the Standard Model expectation, keeping the Higgs production 
rate close to the Standard Model prediction. Larger deviations of A t are forbidden by constraints on the ttH channels. 
For example, the combined CMS analysis of ttH channels reports a signal strength of 2.8+J'g [32], and a specific CMS 
analysis of the ttH —> ttbb channel arrives at a signal strength of 1.2 g [33]. Typical uncertainties around 100% on 
the cross section translate into a 30% uncertainty on the top Yukawa coupling. 

One reason to consider the ttH measurements with care is that their significance hardly adds to an independent 
evidence for this production channel. For example, the combination of Ref. [32] rules out the Standard Model at two 
standard deviations, just slightly less significantly than it establishes the ttH production process. An appropriate 
hypothetical question to ask is if these results would have been published the same way if the signature had been a 
sign of physics beyond the Standard Model instead of a very much expected signal. 

Aside from the large error bars these measurements in the ttH channel have less obvious control of the signal 
kinematics than other Higgs channels; for example, they might or might not include a clear Higgs mass reconstruction, 
which is crucial for the unambiguous interpretation of the rate measurement but poses a well known combinatorics 
problem [43]. Such a global analysis does lead to a valid upper limit on the ttH cross section, but for a lower limit 
we need to assume that ttH production is the only source of relevant events. 

As expected, the individual error bars for A g and A t are around three times as large as the error bar for the 
combination A® M+NP , where the latter is known to better than 20%. The remaining Higgs couplings are again hardly 
affected by the additional parameter A g . The error bar of A 7 is slightly increased because of the enlarged error bar 
on A t . Unlike for A 7 this is a signal for a very strong correlation between A t and A g in the 2-dimensional profile 
likelihood. 

In Fig. [3] we show some relevant 2-dimensional correlations of coupling modifications as obtained for the discussed 
analysis spanning the seven coupling modifications. First, we see that in the A t vs A 7 plane there are four solutions 
corresponding to a sign flip in each of the two couplings. We fix the global sign of all Higgs couplings to Aw > 0. 

As long as we limit our analysis to total rates each individual coupling modification at tree level will show a perfect 
degeneracy between A^ = 0 and A^ = —2. The loop-induced Higgs-gluon coupling is dominated by the top loop, 



Figure 3: Correlations between different coupling modifications for the fit including A 7 as well as A 9 . The 1-dimensional profile 
likelihoods correspond to the results shown as the blue bars in the right panel of Fig. [2] 
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Figure 4: 68% CL error bars on the deviations from all Standard Model couplings of the observed Higgs boson. In addition 
to all couplings predicted by the Standard Model we include a Higgs decay to invisible particles. The results labelled ‘SM exp’ 
assume central values on the Standard Model expectation, but the current data error bars. 


with a small contribution from the bottom quark, so it will not lift this degeneracy. In contrast, the Higgs-photon 
coupling is strongly sensitive to the relative sign of the top and W-contributions. 

The moderate positive correlation in the SM-like solution reflects the fact that an increase of the top Yukawa 
coupling leads to a decrease in the Hj'y coupling and hence has to be compensated by a positive value of A 7 . As 
shown in the central panel the correlation between A t and A g is the strongest correlation in the Higgs couplings 
analysis. It reflects the fact that the gluon fusion Higgs cross section constrains the sum of the two with a slight 
re-weighting from the top mass dependence of the loop-induced Higgs-gluon coupling [44] . We will come back to this 


aspect when discussing the effective theory analysis in Sec. Ill and top mass effects in Sec. [TV] The resulting correlation 
of A 7 and A g first of all features eight solutions, arising from the indirect combination through A t ~ —2, 0. They are 
clearly separated into the two regimes A 7 = —2, 0, while in A g they are merged through the strong correlation with 
A*. For example in the SM-like regime the correlation between the two loop-induced couplings is at a similarly weak 
level as the correlation between A 7 and A t . 

Without showing any detailed results we can also take advantage of the first studies of the Higgs interaction with 
a photon and a Z boson [23; [24] ■ We include a new physics contribution to the loop-induced vertex in the Standard 


Model, in complete analogy to the modifications A 7 and A g in Eqs.(2.1) and (2.2). The corresponding 68% CL allowed 
region on Az 1 bounds A z-y < 0.7 (1.8 at 95% CL), without any visible effect on the rest of studied parameters shown 
in Fig. [2] 


B. Invisible decays 

Higgs decays to invisible particles can only be observed in Higgs production channels with a measurable recoil 
system. Examples are weak boson fusion [35] and ZH production [33], where the more sensitive weak boson fusion 
(WBF) channel might be able to probe invisible branching ratios to 2 — 3% with an ultimate integrated luminosity of 
3000 fb _1 [45]. To date there are ATLAS and CMS analyses available in these two channels E5H5U] . 

In Fig. [4] we show the status for the full set of SM Higgs couplings and a hypothetical Higgs coupling to invisible 
states. Unlike for the other couplings we do not define a coupling deviation A; nv , but directly refer to the invisible 
branching ratio BRi nv . In the Standard Model this invisible branching ratio is generated by the decay H —► ZZ* —> 

It only reaches around 1% and is therefore unlikely to ever be observed at the LHC. The current limit on invisible 
Higgs decays in the full Higgs couplings analysis is around 10%. Obviously, for a dedicated analysis with a more 
constraining model assumption the limits will be stronger [7]. 
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Both relevant production processes responsible for invisible Higgs decay searches are mediated by the ZZH and 
WWH interactions, where A wz are the best measured couplings in the analysis described in Sec. IIA The additional 
searches for invisible decays will not add any new information on the determination on A w,z, so we expect the invisible 
branching ratio to be orthogonal to the other Higgs coupling measurements, i.e. uncorrelated with all other channels. 
A slight correlation of the invisible contribution to the total Higgs width leads to a minor upwards shift of all other 
couplings. 

In Fig.[5]we show the 2-dimensional profile likelihoods for BR; nv versus the Aw and A z appearing in the production 
processes and with A;, dominating the total width. None of them show a significant correlation. In the absence of 
strong correlations with any other model parameters in the Lagrangian of Eq. (2.2), our best fit value and 68% CL 
limits on an invisible Higgs branching ratio BR; nv = 0.16 ^q'^ will hardly depend on the assumptions for example 
made about the loop-induced Higgs couplings. 


C. Theoretical uncertainties 

In this last part of the Higgs couplings analysis we highlight open questions related to the treatment of theoretical 
uncertainties. Unlike experimental uncertainties, the estimate of for example higher-order contributions missing in 
the calculation of an LHC cross section do not offer a frequentist interpretation. It is unclear what kind of likelihood 
distribution of the associated nuisance parameters we need to assume. On the other hand, once we define a likelihood 
distribution for these nuisance parameters the rest of the likelihood analysis is completely defined. 1 

All we can say from a theory perspective is that a certain deviation from the best available cross section or rate 
prediction is in some kind of agreement with the Standard Model or beyond the level where we are willing to consider 
such an interpretation. This problem is independent of the way we determine the uncertainty range on an observable. 
Varying the unphysical factorization and renormalization scales is only one method, and there might be many others. 
If we assume a flat distribution for the theoretical uncertainty to remove any bias between different predictions within 
the allowed range, the RFit scheme m is uniquely defined as the profile likelihood combination with the experimental 
uncertainties. Uncertainties on the parton densities are treated in complete analogy to the theoretical uncertainties 
from unknown higher orders, including any assumption on their correlations. 

Note that a similar problem of choosing a prior for the theoretical uncertainty arises in the Bayesian approach [48] . 
On the one hand the Bayesian approach allows for a choice of priors in general, including the theoretical uncertainty, 
without having to ask for a statistical interpretation. On the other hand, this renders one assumption on the prior as 
ad-hoc as any other. This leaves us with the crucial task to carefully check the prior dependence of our result. 



Figure 5: Correlations between different coupling modifications to SM particles and the invisible branching ratio. The corre¬ 
sponding 1-dimensional profile likelihoods are shown as the blue bars in Fig. [4] 


t Because of this lack of uniqueness in the definition of theoretical uncertainties we advocate for not including them in the experimental 
analyses or (if unavoidable) for factoring them out to allow for a flexible analysis m- 
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No matter what approach we follow, it is important to recognize that we should employ a conservative estimate of 
theoretical uncertainties in addition to a flexible framework which allows us to test different assumptions as efficiently 
as possible [40]. This includes the shape of the likelihood associated with the theoretical error bars as well as the 
size of the theoretical error bars. Before we discuss the modelling of the likelihood including a theoretical uncertainty 
we need to understand an effect which we observe all through Sec. |II A| and |IIB| the expected size of the error bars 
is consistently larger than the observed errors. The reason is linked to the behavior of flat theoretical uncertainties 
once the measurements start developing a pull. In Fig. [6] we show a set of 1-dimensional profile likelihoods. For the 
expected limits, i.e. assuming that all rate measurements agree perfectly with the Standard Model, we clearly see the 
flat central range, induced by the theoretical uncertainties. 

However, once we allow for a statistical distribution of the measurements all 1-dinrensional profile likelihoods lose 
the flat central regions and instead follow the Gaussian shape of the dominant experimental uncertainties. Note that 
this does not have to be the case based on first principles: if all uncertainties were flat in the rates, the resulting 
profile likelihood would keep its box shape, and the errors would be added linearly mm- The central limit theorem 
does not guarantee a Gaussian distribution, because the profile likelihood does not involve a convolution. The 
resulting Gaussians in Fig. [6] instead reflect the fact that theoretical uncertainties are smaller than their experimental 
counter parts, and the Gaussian features of the latter dominate the final distribution once we allow for a spread of 
measurements. The curves in Fig. [6] illustrate the general observation, that actual 1-dimensional error bands with 
their Gaussian behavior are smaller than the expected errors with their flat central range, once we include real data. 

With this observation in mind we show a set of results from our systematic study of the appropriate treatment of 
theoretical uncertainties in Fig. [7] The SFitter standard approach to theoretical uncertainties is based on 

• uncorrelated uncertainties for the production, to account for very different kinematic selections; 

• an allowed cross section range given by the scale dependence of the best available prediction; 

• a flat likelihood distribution of the associated nuisance parameter for the cross sections and the decays. 

These three assumptions we check one-by-one. In the upper left panel of Fig. [7] we show an alternative SFitter 
analysis with fully correlated uncertainties for the production rates, including the error bar from the parton densities. 
Because of the strong correlations between the different production and decay processes expressed in terms of Higgs 
couplings such a shift in the assumed correlations for the theoretical uncertainties could have a significant effect. 
However, we see that for full correlations the size of the error bars is only slightly reduced, in spite of the fact that 
the central values for example for A h shifts by half a standard deviation. The truth lies somewhere in between the 
fully correlated and the fully uncorrelated theoretical uncertainties, where for the upcoming Run II there will be a 
tendency towards less correlation because of the more specific analysis strategies. On the other hand, the difference 



Figure 6: 1-dimensional profile likelihoods for Aw, A z, and A® M+NP . We assume all measurements on the Standard Model 
values with flat theoretical uncertainties (black), observed rates with flat (solid red) and Gaussian (dashed red) theoretical 
uncertainties. 
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Figure 7: 68% CL error bars on the deviations A x from all Standard Model couplings of the observed Higgs boson. First, we 
show the effect of fully correlated theoretical uncertainties on the different production processes (upper left); next, we show the 
results with the modified theoretical uncertainties proposed in Ref. [50] (upper right); then, we illustrate the effect of the N li LO 
calculation of the Higgs cross section [51] (lower left); finally, we illustrate what happens when we simulate the theoretical 
uncertainties with a Gaussian distribution (lower right). The results labelled ‘SM exp’ assume central values on the Standard 
Model expectation, but the current data error bars. 


between fully uncorrelated and fully correlated errors is not worrisome and we stay on the conservative, uncorrelated 
side. 

In the upper right panel of Fig. [7] we show the change in the extracted Higgs couplings when we modify the assumed 
theoretical uncertainties following Ref. m- For a collider energy of 8 TeV and a Higgs mass of 125 GeV the default 
prediction of the gluon fusion production cross section is 

a pp ^H = 19.52 pb ± 7.5%pert ± 7.2% pdf . (2.5) 

Using the modified definition of Ref. [50] it becomes 

<jp P ^H = 22.49 pb ± 13% P ert ± 7.2% p df. (2.6) 

The change in central values as well as in the size of the error bars is hardly observable, in spite of the sizeable change 
in the size of the theoretical uncertainties. This confirms the earlier observation that theoretical uncertainties are 
subleading for the Run I analysis. 

Also related to the theoretical uncertainty we include the recent computation of the N 3 LO corrections to the Higgs 
production rate at the LHC m- The corresponding results are shown in the lower left panel of Fig. [7} Now the cross 
section prediction reads 


Opp^H = 19.95 pb ± 1.5% p ert ± 7.2% p df 


(2.7) 
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with a significantly more optimistic error based on the central scale choice hr,f = m-H / 2. Again, the now strongly 
reduced theoretical uncertainty hardly affects the Run I results. 

Finally, we compare the precision of the Higgs couplings determination with flat theoretical uncertainties with a 
Gaussian nuisance parameter. The main differences between the frequentist RFit treatment and Gaussian theoret¬ 
ical uncertainties are not related to the shape of the final distribution, but to the size of the combined theoretical 
uncertainties. First, combining two flat theoretical uncertainties, for example from unknown higher orders and the 
parton densities, will lead to a linear combination of the two error bars in the frequentist RFit scheme [471 SSI [52]. 
In the Gaussian approach they are added in quadrature. Second, it is not clear with which Gaussian significance we 
should identify the ends of the box-shaped distribution. For example, computing the standard deviation of a flat data 
set stays well below the size of the box. This means that if we compare the range of one standard deviation for the 
RFit scheme with one standard deviation of the Gaussian, the error on the flat distribution appears smaller. 

Per se, it is not clear which of the two effects will dominate in a given fit. In this situation we could choose the flat 
and Gaussian theoretical uncertainties without a clear preference. We stick to the former because we assume that it 
will be the conservative approach once theoretical errors actually affect LHC results with larger data sets. 


III. HIGGS OPERATORS 

Going beyond a measurement of all couplings predicted by the Standard Model we can ask a different question: 
Which consistent Lagrangian describes all LHC measurements best? A standard approach is defined by effective field 
theory [53] . where we categorize a Lagrangian with the appropriate symmetries in terms of the expansion parameter. 
While the results of the previous section can be interpreted in the framework of a non-linear effective Lagrangian 
approach as we have explained, in this section we focus on the linear case. In the linear sigma model we construct a 
SU(2)l x [/(l)v-symmetric Higgs Lagrangian based on the doublet <I> and order it according to the inverse powers of 
the cutoff scale, 1/A [321 IS'IHSl] . The Lagrangian, here restricting to all dimension-6 operators 

C = (31) 
X 

is gauge invariant, but not fully renormalizable or unitary. 

Strictly speaking, in the SM Higgs sector we should separate two sources of dimension-6 operators. Yukawa 
couplings or gauge boson couplings from spontaneous symmetry breaking violate the Appelquist-Carazzone decoupling 
theorem m, which means that the Higgs couplings to photons and gluons are only suppressed by 1/v. New physics 
generally gives rise to dimension-6 operators suppressed by 1/A 2 , leading to Higgs coupling strengths to photons and 
gluons scaling like v/A 2 . This distinction will be reflected in the normalization of the respective operators below. 


A. Dimension-6 operator basis 

Before we present the result of the LHC analysis we need to define our basis of dimension-6 operators. The 
minimum independent set of dimension-6 operators with the SM particle content (including the Higgs boson as an 
SU(2)l doublet) and compatible with the SM gauge symmetries as well as baryon number conservation contains 59 
operators, up to flavor and Hermitian conjugation [51]. To present our choice of operator basis [5], we start by imposing 
C and P invariance and employing for the bosonic sector the classical non-minimal set of dimension-6 operators in 
the HISZ basis [55], with the following operators contributing to the Higgs interactions with gauge bosons: 


Ogg = T'A* G“ v G a ^ 

O ww = <h'W„ v W llv <h 

Obb = &B^B^$ 


Obw = 

O w = (D^W^(D^) 

Ob = {D^B^{D„<f>) 



0 *.2 = \d>‘ ($*$) d „ 

0, PA = (D^) f (£> M $) (<F t $) . 

(3.2) 
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Here the Higgs doublet covariant derivative is D ^4* = (d g + ig'B^/2 + iga a Wf/2) $, the hatted field strengths are 
B gil = ig'B^/2 and W fIl/ = iga a W“ v /2, with the Pauli matrices written as a a . The SU{2) L and U(l)y gauge 
couplings are g and g'. The additional operator £>$,3 = (4A<fi ) 3 is crucial for the structure of the Higgs potential and 
for a theoretical interpretation of the measurement of the Higgs self-coupling, but we can safely omit it for the LHC 
Run I analysis. 

The final choice of structures for our global Higgs analysis follows Ref. [5], relying on operators contributing to 
existing data. We first use the equations of motion (including all necessary fermionic operators }56i omitted in this 
brief introduction) to rotate to a basis where there are not blind directions linked to electroweak precision data. We 
then neglect all operators contributing to the bulk of electroweak precision data at tree level; their coefficients will be 
too constrained to lead to observable deviations in LHC Higgs measurements. After using the remaining equation of 
motion to remove redundancy, we finally neglect all operators that we know will not be constrained by LHC Higgs 
measurements. We are left with the final set of nine operators that parametrize the Higgs interactions at the LHC. 
For the gauge boson interactions they are 


rHVV 


; fcG 


One; + 


fBB 


fww, 


f R 


fxxr 


(3.3) 


The operator 2 appears in the gauge and fermionic Lagrangians, because it leads to a finite renormalization of 
the Higgs field and hence a universal shift of all Higgs couplings to Standard Model fields [9:. This set of dimension -6 
effective operators gives rise to the following Higgs interactions with SM gauge boson pairs, 


-HVV 


= g Hgg HG^G a ^ + g Hll HA gv A^ + g^ A^df'H + HA^Z*” 
+ g$ zz + g {2) zz HZ^ + gf zz HZ g Z» 

+ 9hww (W+ W H + h.c.) + g\ 1 \ vw HW+ V W M + g^ww HW^W 


,( 2 ) 


(3.4) 


where V jlv = d g V v — with V = A,Z,W,G. These effective couplings are related to the coefficients in Eq.(3.3) 

through 


9Hgg — 
QH'yy = 


0‘s fGGV 

8n A 2 

g 2 vs 2 w Ibb + fww 
2A 2 2 


( 1 ) _ 9 , 2 v_ c 2 w fw+s 2 w f B 
9hzz 2A 2 2cl 

( 2 ) _ g 2 V S^fsB + C^, fww 

9hzz ~ 2A 2 2c?„ 


j 2 i 

2A 2 


( 1 ) g 2 v s w (fw - fs) 
9hz 1 ~ 


j 2 i 

2A 2 


2 c,, 

„2 


( 2 ) _ 9 V S w (2sl,fBB — Zc^fww) 

g hz^ — 


2 c,, 


9hww 


g 2 y fw 
2A 2 2 


( 2 ) _ _g?v_ f 

9hww — 2A 2 9ww 


9hzz = m 2 z (V2G F ) 1/2 ^1 - 2 ^ 9hww = m 2 w {V2G F ) 1/2 ^1 - 


(3.5) 


where s w and c w stands for the sine and cosine of the weak mixing angle. 


We finally focus on the huge set of dimension -6 operators contributing to the Higgs interactions with fermion 
pairs [55) . Because of a lack of appropriate observables in the LHC Higgs measurements, from the fermionic operators 
left in the final basis we limit ourselves to the flavor-diagonal Yukawa structures 

O e <S>,33 = ( < fi^)(A3 < i ) e_R ! 3) 0„<6, 33 = (4^4>) (Q3&UR,3) Od<5 >,33 = ( < f’^4>)((33 < i>dfl i 3) , (3.6) 


with $ = j<72<i>*, and where the conventions for the fermion fields are L for the lepton doublet, Q for the quark doublet, 
and f F for the SU(2) F singlet fermions. The corresponding effective Lagrangian for the fermionic interactions reads 


r Hff _ 
'-eS ~ 


f T lJl T 

vA 2 


f-^e<3>,33 + 


fbirth 
vA 2 


Od<S>, 33 


ftmt 

vA 2 


Ou$, 33 + 


f<S>, 2 

A 2 


Od>,2 ■ 


(3.7) 
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As mentioned above, C $,2 affects the Higgs couplings universally. In analogy to the Higgs-gluon coupling we scale the 
fermionic f x by a factor m/v to reflect the chiral nature of the Higgs coupling operator [58]. For the Higgs couplings 
to SM fermions this implies 

C Hff = g f Hf L f R + h.c. with g f = ^1 - ^h ,2 ~ > ( 3 - 8 ) 

where we define the physical masses and fermions in the mass basis (/ = r, b , t). 


B. Rate-based analysis 


As a first step we update the global analysis of dimension-6 operators based on the complete Run I data in the 
SFitter framework. The main difference to the analysis of Ref. [9] is the variable top-Yukawa operator, which can 
now be constrained by ttH production as well as the Higgs production via gluon fusion. 

The contributions of the dimension-6 operators to the production rates and decay widths are calculated using 
MadGraph5 [59] and FeynRules [60]. We check our results with Comphep [61] [62] and VBFNLO [63 . We 
approximately include higher-order corrections through AT-factors computed for the Standard Model processes f521j . 
Similarly, for this rate-based analysis we assume that all detector efficiencies are identical for both the SM Higgs 
processes and the corresponding dimension-6 contributions. The results of this 9-parameter global analysis are shown 
in Fig. [8] and Fig. [9] after performing a statistical analysis as described in Sec. TB For the present case we show the 
multiple degenerate solutions. 


In Fig. [8] we depict a selection of interesting correlations between the dimension-6 operators. In addition to the 
correlations discussed in the previous section, e.g. /qg vs ft shown in Fig. [3] the dinrension-6 operators introduce a 
rich structure of correlations related to the Higgs interactions with electroweak gauge bosons. As long as the analysis 
is only based on rate measurements in the Higgs sector, these correlations are the main difference compared to the 
A-framework. The strongest of these correlations is due to the di-photon channel, as it is measured with the highest 
precision. Therefore, the tree-level contributions from fww and /ss to the Higgs coupling to photon pairs generate 
the strong correlation in the left panel of Fig. [8} see Eq.(3.4). The two, slightly separated, allowed regions at 68% CL 
are due to the interference between the dimension-6 amplitudes and the Standard Model ones. The fact that both 
fww and receive their strongest constraints from the di-photon channel, reflects that their contribution in the 
rate based analysis is very similar to the addition of A 7 in the previous section. While this strong correlation is 
partially broken by their smaller contribution to the other channels in the analysis, we will see in the following section 
that the addition of kinematic distributions will increase the sensitive to fww and f R B stemming from VBF and 
Higgs associate production channels. 

In the central panel of Fig.[8]we show the correlation between fs and fw- The Wilson coefficient fw is much more 
strongly constrained than fs, because of the large contributions of the former to the HVV vertices (V = Z 1 W ± ) 



Figure 8: Correlations between different coefficients f x /A 2 , measured in TeV 2 . The 1-dimensional profile likelihoods corre¬ 
sponding to these results are shown as the blue bars in Fig. [9] 
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Figure 9: Error bars on the coefficients f x / A 2 for the dimension-6 operators defined in Eq.( |3.3[ ) (left panel) and Eq.( |3.7| l (right 
panel). We only include total rate information and show 68% CL as well as 95% CL contours. LTnlike for the other 1-dimensional 
profile likelihoods we keep track of the secondary minima in this one figure. The results labelled ‘SM exp’ assume central values 
on the Standard Model expectation, but the current data error bars. 


while f B only contributes to HZZ with a weak mixing angle suppression. The mild impact of f B will eventually be 
compensated by measurements of H —>■ Zj decays. Moreover, the contributions to HVV also correlate f B to fww 
and /bb, as displayed in the right panel of Fig. [8] 

The universal contribution of /$ j2 to all Higgs couplings strongly correlates this operator with the rest of dimension- 
6 structures, both in the bosonic and in the fermionic sectors. This way, 2 in principle lifts the degeneracy between 
the two allowed regions for fb and / T , which is due to the interference between the SM amplitudes and the higher¬ 
dimensional operators. The actual likelihood values for the two minima are still equivalent though. In the A-framework 
these regions are almost entirely degenerate, allowing us to focus on the SM-like solution in that case. 

Starting with the assumption that to first approximation the rate-based analysis of dimension-6 operators is physi¬ 
cally equivalent to the Higgs coupling analysis described in Sec. [HA] the strong correlations shown in Fig. [8] still pose 
a technical problem. The Higgs coupling modifications A x are by definition well aligned with the experimental mea¬ 
surements, which means that the profile likelihood construction down to 1-dimensional likelihoods is straightforward. 
For example the correlation between f B and fw makes it obvious that a profile likelihood either in f B or fw will 
have to deal with strongly non-Gaussian distributions, including secondary minima. 


With this technical caveat in mind we show in Fig. [ 9 ] the best fit points and the corresponding 1-dimensional 68% 
and 95% CL regions for each effective operator. We follow the procedure described in Sec. |IB[ in this case keeping 
all possible solutions for O t ^,r and Ogg- As we have discussed the strongest constraints apply to fww and f BB . 
Next are fw and f ^ 2 , and finally the weaker constraint f Bl as discussed above. Just like for A t the free value of f t 
enlarges the error bars for foe and splits the allowed parameter range into more or less distinct regions, like those 
shown in Figs. [3] and [9] 

In the right panel we observe the expected secondary solutions for all three ft,b,r- To compare the errors on the 
couplings to fermions we should keep in mind that in Eq.(3.7) the chiral factor is taken out of the definition of the 
operator and its associated scale A. As we can see, at the 68% CL the secondary solutions appear as clear additional 
structures, while at 95% CL the SM-like and secondary solutions barely separate for fb . This allows us to cleanly 
separate SM-like solutions from those with merely switched signs of the Yukawa couplings. Note that the latter 
correspond to a new physics scale A ~ 150 GeV in the presence of a chiral symmetry factor, shedding some doubt on 
the effective theory treatment as a whole. 

From a statistical point of view it is not clear how one would deal with such alternative solutions; in our case we 
show the solutions with flipped signs of the Yukawa couplings in Fig. [9j but will omit them in the 1-dimensional profile 
likelihood for the rest of the present section. In the Markov chain analysis they will be of course still included. We 
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will revisit this issue in Sec. IV for the case of the top Yukawa coupling. 


C. Kinematic distributions 


Based exclusively on total event rates, the results from the previous section do not take full advantage of the available 
information. In Eqs.(3.4) and (3.8 1 we observe that 0^,2, O t , and O t merely modify the SM coupling strengths, 
but the other dimension-6 operators do generate new Lorentz structures. These anomalous Lorentz structures are 
best visible in Higgs production rather than decays, because the momentum flow is not limited by the Higgs mass. 
Their study is indeed one of the most interesting aspects of our effective field theory analysis. 

To establish a framework for an implementation of kinematic distributions into the Higgs operator analysis we 
first focus on VH production and weak boson fusion. Adding kinematics to our global analysis faces a considerable 
challenge, because we are limited to fully documented distributions. When multi-variate analysis techniques are 
applied, the documented distributions are usually not optimized. However for two test cases we will show how we can 
consistently combine information from rates with kinematic distributions without weakening the analysis. 

Finally, the effective Higgs Lagrangian does not define a UV-complete theory if we only include dimension-6 op¬ 
erators. The cutoff scale A is encoded in the ansatz, and at least for a weakly interacting theory the experimental 
sensitivity offers a consistency test. There exist several ways to define a model which we can consistently compare to 
data: 


1. take the alternative model at face value, unless a prediction actually violates unitarity. This approach maximizes 
the distinguishing power of the measurement, but it only rules out the ultraviolet completion with the least 
SM-like behavior. 

2. attach momentum-dependent form factors to soften the ultraviolet behavior. The main problem is that after 
going through a lot of trouble of defining an effective field theory hypothesis, we spoil it by introducing ad-hoc 
non-local interactions in the position-space Lagrangian. 

3. only use data in phase space regions which are not sensitive to the ultraviolet completions, for example requiring 
Pt < 100 GeV for the tagging jets in weak boson fusion [64] or an upper bound on p\ in VH production. The 
obvious disadvantage of this approach is that we lose experimental information and produce worse bounds, as 
we will see. 


In this analysis we will attempt to include as much of the kinematic information as possible, but carefully check how 
much of the distinguishing power comes from phase space regions not obviously consistently described by the effective 
field theory. 


We start with VH —> V ( bb ) production. To be maximally sensitive to Owwi Obbi and Ob requires a kinematic 
variable with large flow through the production vertex. A key candidate is the transverse momentum distribution 
in the hard process m of the cut-based experimental analysis, which serves as a check of the measured rate, while 
the measurement itself used in Sec. [II] relies on a multi-variate analysis. The experimental search requires two jets 
with medium and tight &-tags and defines three categories with 0, 1 and 2 leptons that receive contributions from 
HW and HZ productions, therefore being sensitive to different contributions of the higher-dimensional operators. 
We show these three distributions in Fig. [TO] with a selection of dinrension-6 anomalous contributions. While the 
background rapidly decreases at large transverse momenta, the main effect of the dimension-6 operators is conversely 
an enhancement at high momenta, being most prominent in the last bin of the distributions. 

To combine the different lepton multiplicities we use our FeynRules ESI implementation of the dimension-6 
operators to generate the distributions with MadGraphS [55], Pythia [55], and PGS4 [BB], the latter checked with 
Delphes Ej. We use the SM Higgs expectations to calibrate our setup to the distributions and rates shown in 
Ref. [21] . We parametrize the kinematic distributions as a function of the dimension-6 operators and we use the SM 
background expectations and the number of measured events per bin in Ref. m ■ 
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Figure 10: Upper and left lower panels: Pt distributions from VH production |2T|. We show the SM Higgs-plus-background 
expectation (black solid), the number of observed events, the SM Higgs hypothesis (red solid), and the expectation from one 
dimension-6 operator (red dashed). Lower right: A <j>jj distribution in Ref. [34]. We display the number of observed events, the 
SM-Higgs hypothesis (red), and the expectation from adding a set of dimension-6 operators (dashed red and dotted blue). All 
plots include 20.3 fb -1 at 8 TeV. In the figure we neglect the effect of higher-dimensional operators on the branching ratios. 


As mentioned above, the pj' distributions shown in Fig. 10 correspond to a cut-based ATLAS analysis [21]. When 


adding this information to the global Run I analysis consistently we need to be careful: first, we cannot use the same 
information twice. This means we could remove the corresponding total rates from the analysis and instead include 
the binned distributions. However, the cut-based analysis is weaker than the multi-variate analysis and they do not 
give the same measured central values. This would render any estimate of the additional power of the kinematic 
information impossible. Instead, we keep the multi-variate rate information and add the kinematic information 
through a set of asymmetries based on the bin content of Fig. [ToJ 


A, = 


bin .; +1 - bin, 
binj+i + bin.; 


(3.9) 


which for each leptonic channel defines three or four additional measurements. 


Our second test case is the azimuthal angle correlation in weak boson fusion production with H —> 77 [55]. Because 
the measurement of A <pjj does not require the reconstruction of any reference frame, its uncertainties are reduced. 
Unfortunately, the corresponding distributions are not shown in the most prominent weak boson fusion channels with 
decays H —> W + W~ and H — >• t + t~ . An unfolded distribution is in contrast available for the decay H —► 77 fTTT| . 
However, due to the lack of cuts on rrijj and A r]jj, the weak boson fusion mode accounts for less than 35% of all 
signal events, diluting consequently the promising power of the A <fjj variable in this production mechanism. 

In the present absence of a better alternative we include the above channels in our SFitter analysis. To simulate 
SM Higgs production in weak boson fusion and the VH channel we rely on the same selection of tools we have used 
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for the pj- implementation. To validate our calculations, we compare our SM simulations to the ATLAS result, most 
notably the plots available in HEPDATA [S5J. Once our setup is tested we simulate the effect of dinrension-6 operators 
on the weak boson fusion and VH distributions. The main contribution from Higgs production in gluon fusion is only 
affected by Oqg, O^^, and Ob,t, none of which change the Lorentz structure of the hard process. We can then use the 
central estimate by ATLAS, properly reweighted by the introduced shift of the relevant operators. We use a similar 
reweighting to simulate the effects of the effective operators in the di-photon decay, as none of the operators generate 


a non-SM Lorentz structure for this vertex either. In the lower right panel of Fig. 10 we show the A <fijj distribution 


with a selection of dimension-6 contributions. It turns out that all Oww, Obb, Off, and Ob peak at 0 or at n [SB] , 
To add A.(f>jj to the global analysis we again keep the measured total rates used in the previous analyses and 
construct three additional asymmetries [68j . 


Ai = 


A? — - 


A y, = 


g (A fa < f) + oj Aiftj-j >f)~ j < N>n_ < -f) 

ct (A <t>jj < f) + ff(A (t>jj > f ) + a (§ < < f) 

g (A <t>jj >f)- g(A4> 3 -j < f) 

> x) + ct ( A< Ah < f) 

a(Acj }jj > 5 J) + a(^<A(j> jj < 5 J) ' 


(3.10) 


The first asymmetry is tailored to discriminate different production modes and CP structures |34l I68j . The second 
asymmetry enhances the sensitivity to Oww, Obb , Ow, and ®b in weak boson fusion, all of which generate non-zero 
values for A 2 . Finally, the third asymmetry is orthogonal to the other two, to not exclude any information. 


D. Full dimension-6 analysis 

The final step in our higher-dimensional SFitter analysis is to add these test distributions to the coupling infor¬ 


mation used in Sec. IIIB As we have discussed, we use the experimental information shown in Fig. 10 [21, 34] in 


terms of the asymmetries defined in Eqs.(3.9l and (3.10). 


The main technical problem of the purely rate-based analysis of dimension-6 operators are the correlations which 
make it hard to extract 1-dimensional profile likelihoods and error bars for the individual f x / A 2 , as illustrated in 
Fig-! We show the effect of kinematic distributions on some critical 2-dimensional profile likelihoods in Fig. |TT] 
In the top row we show the results after including the A <f>jj distribution only. Compared to the Fig. [8] we see very 
small improvement, except for a slight reduction of the secondary structure in the Obb vs Ob and Ow vs Ob 
correlations. However, this reduced impact should not be taken as a statement about the distinguishing power of the 
A <f)jj distribution; it is really linked to the lack of publicly available information on this distribution, as discussed 
above. 

we show the impact of also adding the full p}f information in VH production. It 


In the second row of Fig. 
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significantly improves the situation with secondary solutions, largely removing the correlated structure for example in 
Ow vs Ob- There still exists a weak secondary minimum for example in the Obb vs Ob correlation, but because of its 
relative weakness it will allow us to derive a more straightforward 1-dimensional profile likelihood and an associated 
68% CL error bar for example on /bb/A 2 . The reduction of the allowed space is notorious in the three corresponding 
panels in Fig. 11 Actually, the left panel shows that after including the p\- and Afyj distributions the fit becomes 


more sensitive to Oww and Obb through their individual contributions to the HVV couplings mediating associated 
VH production and weak boson fusion. 


Finally, we need to check the consistency of the effective theory approach m ■ Based on Run I data our analysis 
typically probes \f x /A 2 \ ~ 10/TeV 2 . A hypothetical setting of f x to unity would correspond to new physics scales 
around A ~ 300 GeV. According to Fig. 10 the highest momentum bin of the p x T distribution starts from pj, = 200 GeV 


and includes all events above this value. A conservative approach would be to exclude this last bin, and thus the last 


asymmetries defined in Eq. (3.9), from the kinematic analysis. In the bottom row of Fig. 11 we show the corresponding 
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Figure 11: Correlations between different coefficients f x / A 2 (measured in TeV~ 2 ) after including kinematic distributions. In 
the top row we add the A (f>jj distribution; in the second row we also include p\ from VH production; in the bottom row we 
then remove the highest bin associated with large momentum flow through the dimension-6 vertex. The 1-dimensional profile 
likelihoods of the second row correspond to the results shown as the blue bars in Fig. |12| 


2-dimensional correlations from this analysis. A comparison to the first row shows that almost the entire additional 
information of the distribution is encoded in the last bin. In the remainder of the discussion we will not follow this 
conservative approach, so it should be noted that the full SFitter analysis of the higher-dimensional operators has 
to be taken with a grain of salt. On the other hand, the analysis including kinematic distributions is mostly meant 
to be a proof of principle, and the consistency of the Higgs effective theory will clearly improve with Run II data. 

With this in mind we show the best fit points and the corresponding 1-dimensional 68% CL error bars including 
kinematic distributions in Fig. [l2j In contrast to Fig. [9] we do not show secondary solutions for the signs of the 
Yukawa-like couplings and for Oqg- We also limit ourselves to 68% CL contours. We see that Ob and Ow are the 
operators most affected by the addition of kinematic distributions, closely followed by Oww and Obb- Typical energy 
scales probed by Run I data are 300 GeV to 500 GeV if order one Wilson coefficients are assumed, with less significant 
constraints in the fermion sector. All coefficients are in agreement with zero, and the one to two sigma deviations 
are hard to map onto individual measurements. Including all available kinematic information visibly stabilizes the 
constraint on fs and moves every single best-fit point closer to the Standard Model prediction. 


IV. FUTURE: OFF-SHELL MEASUREMENTS 

ATLAS and CMS recently published a study on the contribution of Higgs exchange to ZZ production at invariant 
masses well above the Higgs pole mzz ~ m# DU E3 E5J. Given the small Higgs width, such a measurement 
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Figure 12: 


CL error bars on the Wilson coefficients f x / A 2 for the dimension-6 operators defined in Eq.pOJ (left panel) 


and Eq.(3.71 (right panel). In addition to total rate information we also include kinematic distributions and only show 68% 
CL contours. For the Yukawa couplings as well as for Ogg we limit ourselves to the SM-like solution for this representation. 


would normally only show a very moderate dependence on the Higgs mass. However, the kinematic structure of 
this particular channel turns it into a sensitive measurement. Approximately 0(15%) of the rate mediated by the 
s-channel Higgs exchange lies in the off-shell regime, mu > 130 GeV. In addition the leading effect arises from the 
signal interference with the continuum background. Some representative Feynman diagrams to this process are 



V 





Note that A g in this representation shows a non-trivial momentum dependence, limiting the model-independent 
features of the width measurement m■ If the Higgs propagator in the interference is probed far above the mass shell, 
it behaves like 1/s. On-shcll and off-shell Higgs rates then scale like 


on-shell 


OC 


gU m H) g 2 f (m H ) 

T h 


vs CT off-shen ^ ^ g 2 ( TO4<! ) gj ( m4€ ) . 


(4.1) 


where gi ( gf ) refer to the Higgs couplings involved in the production (decay) for the present channel. Eventually, 
we will remove the assumptions about the Higgs width described in Sec. |I B| from the SFitter setup and instead 
determine the total width from the combination of off-shell and on-shell measurements. The Lagrangian of the 
underlying hypothesis reads 

£ = £sm + Aw gmwH + A z m z H ^ A/ — -H (Jr/l + h.c.) 

ZCW T,b,t V 

H H 

+ A „Fg — GnvG^ + A jFa — + invisible decays + unobservable decays . (4.2) 

v v 

The distinction between the two terms linked to decays to non-SM states is that ‘invisible decays’ are reconstructable 
using missing transverse momentum, while ‘unobservable decays’ are for some other reason not observable at the 
LHC, for example because of overwhelming jet backgrounds [5]. Not accounting for such unobservable decays would 
lead to shifts of all A^ as compared to the analysis including these decays. 


Before we allow for a fully unconstrained Higgs width through unobservable decay channels we combine on-shell 
and off-shell analysis to probe the energy dependence of the operators involved EP- On the Higgs production side, the 
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m 4 i [GeV] 


L=20.3 (8 TeV) fb ' ATLAS 



m 4 | [GeV] 


Figure 13: mu distribution for the CMS (left) and ATLAS (right) analyses. The qq —> ZZ background and the data points 
are obtained from Refs. [35l l36| . The remaining curves are generated following the parametrization of Eq.(4.31. 


dimension-6 operators entering the off-sliell measurements are listed in Eqs.(3.3I and (3.7), namely Oqg , Oa >,2 and 
They can be described by the two parameters A 9 and A t or equivalently gHgg and gj (with a marginal contribution 
from the bottom loop). The difference between the two are top mass effects in the kinematic structure [44] . In the 
decay the climension-6 operators in Eq.(3.4) lead to additional operator structures, namely Z^Z^d^H and HZ I , IJ Z I ‘ 1 ' 


in Eq.(3.4). None of them affects the longitudinal Z-polarization m, so they lead to similar mu kinematics as the 


(3) 

SM-operator HZ^Z 9 ’. We parametrize this Higgs decay only accounting for A z or g y HZZ - Properly accounting for 
the continuum background we can write the gluon fusion component to the signal as 


■Mgg^tZZ — (1 + A z ) [(1 + A t)M-t + A 9 A4 9 ] + M. c 


da 

dmu 


= (1 + A z ) 


(1 + A t ) 


da tc 

dmu 


^ da gc 
9 dmu 


+ (1 + A Z ) 


2 < 1 + a '» 2 ^ + < 1 + A -> a * 


dvtg a 2 da gg 
9 dm. 


dmu 


At 


+ 


da c 
dmu 


(4.3) 


We illustrate the top mass effects in Fig. [13] The background qq —> ZZ and the data points are taken from the 
experimental publications [35] [36]. The gluon-initiated component is generated with Mcfm m following Eq.(|4.3|). 


QCD corrections to the gluon-induced component are accounted via a global AT-factor fm 173] . We follow the ATLAS 
and CMS cut-flow and find full agreement with both studies. 



Figure 14: Correlation between the coupling modifications At and A g without (left) and with (center) off-shell Higgs measure¬ 
ments. In the right panel we show the 1-dimensional profile likelihood for At with and without off-shell measurements. 
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Figure 15: Correlation between the total Higgs width and the coupling modifications A w (left) and 1-dimensional profile 
likelihood of the total Higgs width (right). Both, on-shell and off-shell Higgs rates are included. 


Following Eq.(4.3), terms linear in (1 + A*) are sensitive to the sign of the top Yukawa coupling [44] : in the Standard 
Model the off-shell interference is destructive, while a sign change in the top Yukawa coupling increases the combined 
rate significantly. In addition, we see the kinematic difference from the missing top mass threshold and the missing 
logarithmic top mass dependence for the dimension-6 operator. Similar top mass effects can be observed in gluon 
fusion Higgs production with hard jets (21 jJJJ] and in the gluon-induced contribution to VH production [75]. Both 
should eventually be included in the Higgs couplings analysis. 

One way to exploit this feature in our Higgs couplings determination is to include the usual coupling modifications 
and an invisible branching ratio, but no unobservable width. In this case we probe the momentum dependence of the 


effective Higgs-gluon coupling, linked to its top mass dependence |44| . In Fig. 14 we present the resulting correlation 
between A t and A g , finding a significant improvement from the off-shell rate measurement. The SM-like and flipped- 
sign solutions clearly separate. In the right panel of Fig. [14] we observe a slight preference towards a negative top 
Yukawa coupling. It arises from a small excess of events in the off-sliell CMS data. ATLAS sees the opposite trend, 
but with considerably fewer events. 


Finally, we can allow for unobserved decays according to Eq.(4.2). This corresponds to a free total Higgs width 


under the condition that the minimum width is given by the sum of the observed partial widths. In Fig. 1 15] we present 
first results from a global SFitter analysis. In the left panel we show the correlation between a typical coupling Aw 
and the total Higgs width. For r^/Tf^ 1 the Higgs production and decay rates scale like g^/Tn- Indeed, we see 
that the positive correlation extends to r^/T™ ~ 30 ~ 2.3 4 . The corresponding value of Aw is then 1.3, just as 
expected. 

In the right panels we see that the upper bound is approximately Th < 9.3 T™ at 68% CL. While our width 
constraint was obtained in the EFT context, our bound is still competitive to other analysis that account only to 
SM-like interactions m- The key ingredient here is the analysis of the whole m^i distribution profile that probes the 
possible new physics at different energy scales. A similar study via on-shell signal strengths would clearly not be as 
sensitive. 


V. PRESENT AND FUTURE 

In this paper we have presented a final analysis of the Run I Higgs measurements. For the first time, we directly 
compare the direct Higgs coupling analysis with the effective Lagrangian approach. For the coupling analysis we 
included independent variations of the Higgs-top and Higgs-gluon coupling, as well as an invisible Higgs branching ratio 
and eventually Higgs decays through unobserved channels. While for the former the current results for ttH production 
should be taken with a grain of salt, we found a stable combined measurement of BRj nv = 0.16l[lo ii • Including off-shell 
Higgs measurements the full couplings fit gave an upper limit on the total Higgs width of T# < 9.3 rf^ at 68% CL. 
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Theoretical uncertainties and their effect on the individual error bars of more fundamental parameters are starting to 
become a crucial issue already by the end of Run I. We tested different assumptions on the correlation of theoretical 
uncertainties, on their size, and on their statistical treatment. As in the absence of a clear definition based on 
statistical principles we at least ensured that our analysis was conservative. Based on our findings the issue of 
theoretical uncertainties and their statistical treatment needs to be carefully considered for the upcoming Run II 
data. 

To include kinematic distributions the traditional coupling analysis has to be expanded. We consider an effective 
field theory approach the most natural and the most promising expansion. As long as we only include Higgs rate 
information the Run I analysis in terms of Higgs coupling modifications A x and in terms of higher-dimensional 
operators /^/A 2 are essentially equivalent. Once we included measurements from the pure gauge sector this will 
change. Technically, the rate-based analysis in terms of an effective field theory is more challenging, because the 
underlying parameters are less directly linked to measurements than for the Higgs coupling modifications. 

We added two sample distributions from VH production and from weak boson fusion to our analysis. We found 
a significant stabilization of the higher-dimensional analysis, the secondary structures were more easily identified, 
and all individual error bars on the Wilson coefficients were visibly reduced. We end the discussion with a word of 
warning: given the current precision in the dinrension-6 analysis it is not guaranteed that the effective Lagrangian 
expansion is within its range of validity for all considered observables. It should then be seen as a motivated and 
useful parametrization of Higgs interactions, while we wait for the increase of precision in the future Run II of the 
LHC. 
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Appendix A: Numerical results 


In this appendix we give the limits presented in the figures throughout the paper as numbers. 



5 parameter analysis 

6 parameter analysis 

Best fit 

68% CL intervals 

95% CL intervals 

Best fit 

68% CL intervals 

95% CL intervals 

A w 

-0.113 

(-0.23,0.02) 

(-0.335,0.145) 

-0.185 

(-0.35,-0.05) 

(-0.465,0.095) 

Az 

0.0563 

(-0.07,0.18) 

(-0.20,0.295) 

0.041 

(-0.08,0.175) 

(-0.22,0.285) 

A t 

-0.271 

(-0.38,-0.08) 

(-0.495,0.09) 

-0.271 

(-0.41,-0.08) 

(-0.53,0.09) 

A„ 

-0.291 

(-0.535,0.035) 

(-0.785,0.36) 

-0.304 

(-0.57,0.015) 

(-0.83,0.355) 

A T 

-0.0987 

(-0.265,0.095) 

(-0.4,0.26) 

-0.0826 

(-0.265,0.095) 

(-0.405,0.29) 

Ay 

— 

— 

— 

0.129 

(0.015,0.29) 

(-0.11,0.42) 

a sm+np 

— 

— 

— 

-0.033 

(-0.175,0.13) 

(-0.305,0.31) 


(-21nL) min = 69.2, (-21nL) SM = 72.1 

(-2 In L) min = 68.1, (—21nL)g M = 72.1 


Table I: Best fit values, 68% CL and 95% CL allowed ranges for the results of the Higgs analysis with 5 free couplings (blue 
bars in Fig. [T]) and the results of the analysis including in addition A 7 (blue bars in the left panel in Fig. |2|. 



7 parameter analysis 

8 parameter analysis 

Best fit 

68% CL intervals 

95% CL intervals 

Best fit 

68% CL intervals 

95% CL intervals 

A w 

-0.160 

(-0.335, -0.05) 

(-0.46,0.085) 

-0.0867 

(-0.265,0.025) 

(-0.38,0.155) 

Az 

0.0559 

(-0.07,0.195) 

(-0.205,0.305) 

0.158 

(0.01,0.28) 

(-0.125,0.405) 

A t 

0.159 

(-0.2,0.46) 

(-0.585,0.75) 

0.188 

(-0.13,0.57) 

(-0.505,0.845) 

A b 

-0.265 

(-0.565, -0.01) 

(-0.82,0.295) 

-0.193 

(-0.5,0.06) 

(-0.77,0.375) 

A T 

-0.0492 

(-0.25,0.095) 

(-0.395,0.28) 

0.0417 

(-0.17,0.185) 

(-0.33,0.375) 

Ay 

0.226 

(0.09,0.40) 

(-0.065,0.555) 

0.248 

(0.1,0.435) 

(-0.055,0.595) 

A g 

-0.479 

(-0.83,-0.125) 

(-1,0.37) 

-0.430 

(-0.855,-0.13) 

(-1,0.385) 

BRinv 

— 

— 

— 

0.157 

(0.048,0.226) 

(0., 0.306) 

a sm+np 

£A 7 

-0.0191 

(-0.17,0.125) 

(-0.295,0.285) 

0.0892 

(-0.09,0.22) 

(-0.22,0.395) 

a sm+np 

^9 

0.230 

(-0.4,0.115) 

(-0.51,0.35) 

-0.163 

(-0.335,0.04) 

(-0.45,0.115) 


(—2InL) min = 66.4, (-21nL) SM = 72.1 

(—2InL) min = 63.4, (-21nL) SM = 72.1 


Table II: Best fit values, 68% CL and 95% CL allowed ranges for the results of the analysis with 7 free couplings (blue bars in 
the right panel in Fig. [2]) and the results of the analysis including in addition BRs nv (blue bars in Fig. |4|. 



Correlated theoretical uncertainties 

N 3 LO gluon fusion prediction 

Best fit 

68% CL intervals 

95% CL intervals 

Best fit 

68% CL intervals 

95% CL intervals 

A w 

-0.137 

(-0.26,-0.005) 

(-0.39,0.115) 

-0.162 

(-0.315,-0.04) 

(-0.44,0.09) 

Az 

0.0814 

(-0.02,0.225) 

(-0.14,0.335) 

0.0704 

(-0.07,0.19) 

(-0.205,0.295) 

At 

0.113 

(-0.07,0.525) 

(-0.315,0.805) 

0.165 

(-0.195,0.475) 

(-0.565,0.76) 

A;, 

-0.317 

(-0.51,-0.01) 

(-0.725,0.28) 

-0.271 

(-0.55,0.005) 

(-0.8,0.295) 

a t 

-0.043 

(-0.195,0.105) 

(-0.32,0.27) 

-0.0939 

(-0.25,0.085) 

(-0.395,0.26) 

Ay 

0.213 

(0.105,0.375) 

(-0.015,0.52) 

0.241 

(0.09,0.395) 

(-0.06,0.54) 

A a 

-0.386 

(-0.85,-0.215) 

(—1., 0.08) 

-0.508 

(-0.82,-0.11) 

(—1., 0.39) 

A.SM+NP 

aa 7 

0.0088 

(-0.115,0.15) 

(-0.23,0.29) 

-0.00786 

(-0.16,0.12) 

(-0.28,0.28) 

^SM+NP 

-0.195 

(-0.355,-0.055) 

(-0.45,0.11) 

-0.248 

(0.37,-0.095) 

(-0.49,0.045) 


(—2 In L) min = 95.6, (-21nL) SM = 105.3 

(—2 In L) min = 71.8, (-21nL) SM = 77.3 


Table III: Best fit values, 68% CL and 95% CL allowed ranges for the results of the analysis with 7 free couplings assuming 
correlated theoretical uncertainties (light blue bars in the upper-left panel in Fig. [7| and the results of the analysis including 
N 3 LO corrections to the Higgs gluon fusion production rate at the LHC [5Tj (light blue bars in lower left panel in Fig. [t]) . 
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Passarino gluon fusion prediction 

Gaussian theoretical errors 

Best fit 

68% CL intervals 

95% CL intervals 

Best fit 

68% CL intervals 

95% CL intervals 

A w 

-0.174 

(-0.355, -0.06) 

(-0.47,0.09) 

-0.118 

(-0.295,0.01) 

(-0.425,0.14) 

Az 

0.0653 

(-0.08,0.19) 

(-0.21,0.305) 

-0.03848 

(-0.16,0.105) 

(-0.3,0.21) 

A t 

0.139 

(-0.205,0.445) 

(-0.585,0.7) 

0.188 

(-0.235,0.535) 

(-0.835,0.735) 

A& 

-0.291 

(-0.58,-0.015) 

(-0.845,0.31) 

-0.151 

(-0.48,0.18) 

(-0.78,0.585) 

a t 

-0.0674 

(-0.265,0.09) 

(-0.415,0.295) 

-0.0700 

(-0.255,0.1) 

(-0.405,0.305) 

a 7 

0.231 

(0.085,0.4) 

(-0.065,0.56) 

0.150 

(0., 0.315) 

(-0.24,0.45) 

a 9 

-0.552 

(-0.895,-0.205) 

(-1.. 0.305) 

-0.372 

(-0.755,-0.16) 

(-0.96,0.97) 

a sm+np 

^ 7 

-0.0266 

(-0.19,0.115) 

(-0.305,0.305) 

-0.0497 

(-0.215,0.09) 

(-0.36,0.25) 

^SM+NP 

-0.309 

(-0.455,-0.195) 

(-0.55,-0.045) 

-0.117 

(-0.28,0.05) 

(-0.415,0.255) 


(—2 In L) m in = 61.7, (-21nL) SM = 71.0 

(—21nL) m i n = 80.0, (—21nL)sM — 83.0 


Table IV: Best fit values, 68% CL and 95% CL allowed ranges for the results of the analysis with 7 free couplings using the 
gluon fusion prediction from Ref. [50] (light blue bars in the upper-right panel in Fig. [t]) and the results of the analysis assuming 
Gaussian theoretical uncertainties (light blue bars in lower-right panel in Fig. [tJ) . 



Rate-based analysis 

Analysis with kinematic distributions 


Best fit 

68% CL intervals 

95% CL intervals 

Best fit 

68% CL intervals 

95% CL intervals 

fact A 2 (TeV 2 ) 

-19.6 

(-24,16.8) 

(-27.9,-12.9) 

-24.39 

(-27,-18.9) 

(-30,-14.1) 


-4.79 

(-7.2, 0.6) 

(-9.3,9.3) 

-3.75 

(-7.2,-1.2) 

(-9.6, 9.6) 


4.42 

(1.2, 7.2) 

(13.5,28.8) 

3.94 

(0.6, 7.8) 

(15.3,30) 


19.35 

(16.8,25.1) 


24.18 

(19.2,27) 


fww/ A 2 (TeV 2 ) 

4.36 

(-0.4, 9.5) 

(-3.7,13.7) 

0.296 

(-2.95,2.9), (9.65,9.8) 

(-4.3,4.4), (15.3,30) 

/sb/A 2 (TeV 2 ) 

-4.72 

(-8.6, 2.8) 

(-13.4,6.1) 

-0.518 

(-2.45,4.3) 

(-10,7) 

U,2/A 2 (TeV" 2 ) 

6.15 

(0,12.5) 

(-5,18) 

1.03 

(-4.75,6) 

(-8,10), (12.25,19.75) 

fwl A 2 (TeV 2 ) 

-6.38 

(-12.4,7) 

(-13,22.5) 

1.12 

(-2.25,3.75) 

(-5.75,7.25) 

fa/ A 2 (TeV 2 ) 

-29.04 

(-45.2,-5.6) 

(-55.4,13) 

-4.16 

(-12.2,9.4) 

(-45.2,14.8) 

hi A 2 (TeV 2 ) 

4.63 

(-2.1, 7.8) 

(-9.6,10.8) 

0.83 

(-6,6.6) 

(-13.2,19.2) 


33.7 

(24.9,47.3) 

(13.2,61.2) 

46.4 

(33,56.4) 

(24,67.8) 

/r/A 2 (TeV 2 ) 

-3.94 

(-9,2) 

(-12,4.5) 

-1.88 

(-4,3.5) 

(-9,7) 


41.7 

(35,51) 

(24,56.5) 

43.9 

(39,53.5) 

(23.5,60) 

/t/A 2 (TeV" 2 ) 

-5.03 

(-12.8,1.6) 

(-21.3,10.9) 

-4.97 

(-12.9,1.8) 

(-18.5,10.2) 


43.4 

(34,54.4) 

(24.9,63.4) 

53.37 

(41,59.9) 

(24.2,66.9) 


(-2 In 

L) min = 66.7, H 

In T)sm = 72.1 

(—21nL) min = 88.4, (-1 

kiL)sM =91.9 


Table V: Best fit values, 68% CL and 95% CL allowed ranges for the results of the effective Lagrangian rate-based analysis 
(blue bars in Fig. [9]) and the results when including in addition kinematic distributions (blue bars in Fig. |12||. 
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